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ABSTRACT 

Ram pressure stripping of the hot gas that surrounds normal galaxies as they fall 
into groups and clusters (also referred to as 'strangulation' or 'starvation') is generally 
thought to shut down star formation on a time scale of a few Gyr. However, it has 
recently been suggested, on the basis of X-ray-optical scaling relations of galaxies in 
the field and the group/cluster environment, that confinement pressure by the intra- 
clustcr medium can actually lead to an increase in the mass of hot gas surrounding 
these galaxies. We investigate the competition between pressure confinement and ram 
pressure stripping for satellite galaxies in orbit about galaxy groups and clusters us- 
ing simple analytic models and detailed cosmological hydrodynamic simulations. It is 
found that, independent of host mass, ram pressure is generally dominant over con- 
finement pressure — only ~ 16% of galaxies find themselves in the reverse situation. 
Furthermore, these galaxies have, on average, less hot gas than ram-pressure domi- 
nated ones, contrary to simple expectations. This is explained by the fact that the 
small number of galaxies which arc confinement dominated are typically at first or 
second apocentre and have therefore already been maximally affected by ram pres- 
sure stripping around first pericentre. Our results are shown to be insensitive to host 
halo mass; we argue that the same is true for uncertain sub-grid processes, such as 
feedback. 

Key words: galaxies: clusters: general — galaxies: evolution — galaxies: haloes — 
galaxies: interactions — galaxies: intergalactic medium — galaxies: ISM 



1 INTRODUCTION 

In the current paradigm of galaxy formation, haloes of hot 
gas are predicted to be a common feature around galaxies 
with masses similar to (or larger than) the Milky Way 
(|White fc Reesll 19781 : IWhite fc Frenklll99l] h With a temper- 
ature of ~ 10 6 K, they are diffuse sources of soft X-ray emis- 
sion. There is now a large number of observational detections 
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Owen fc Warwickl |2009| ; lAnderson fc Bregmanl 1201 ll : 
Li fc Wand 120121 '). Observed at first only in galaxies in 
the fi eld and poor environments l|Forman. Jones fc Tucken 
Il985h . improved observing facilities like the Chan- 
dra and XMM-Newton telescopes have established the 
presence of hot gas haloes also in galaxies in groups 
dJeltema. Binder fc Mulchaev! 120081 ') and even clusters, 
where Sun et al. (|2007[ ) found them to be as common as 
> 60% in Lk s > "2L* galaxies. 

Studying these X-ray coronae promises to enhance 
our understanding of galaxy formation and evolution, as 
one central predict ion from the theoretical models of 
IWhite fc Reel i|l97Sf ) and lWhite fc Frenkl (|l99lh is that the 
hot gas cools and replenishes the cold gas reservoir which is 
responsible for fuelling star formation. One complication, 
however, is that not only gas which is cooling and thus 
inflowing into the galaxy centre can emit X-rays, but so 
can outflowing hot gas driven by e.g. supernovae or active 
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galactic nuclei (AGN). There is as yet no clear consensus on 
which of these two processes dominates the X-ray emission: 
while it has long been assumed to be the latter (s ee, e.g., 
iRead fc Ponmanlll998l; iMathews fc Brighentilll99g ) , recen t 



work by Crain et al. |2010al . see also Crain et al. l2010bh . 
using hydrodynamical simulations which reproduce the ob- 
served X-ray-optical scaling relations of normal disc galax- 
ies, indicates that the bulk of the X-ray emission comes from 
the cooling of the hot, quasi-hydrostatic corona. Direct con- 
firmation of the nature of the X-ray emitting gas will be 
provided in the future by deep, high-resolution spectra from 
observatories such as IXO or NeXT/ASTRO-H. 

Of course, galaxies are not just shaped by internal pro- 
cesses like gas cooling and outflows, they are also influenced 
by their local environment. For example, the tidal stress 
induced by a group or cluster on an infalling galaxy can 
lead to stripping or even total disruption which may ac- 
count, at least partially, for the well-known morphology- 
density relation of galaxies in dense r environments be - 
ing preferentially of early type (e.g., iMoore et al.l [l99^ ). 
A second important influence is the interaction of the 
intragroup/-cluster medium (ICM) with the galactic gas: 
as the galaxy moves through the ICM, its gas experi- 
ences a drag or ram pressure, which can lead to strip- 
ping o f both the cold, central gas di sc (e.g.. iGunn fc GottJ 
1 19721 : lAbadi. Moore fc Boweil Il999ll and the hot ex- 
tended gaseous halo (e.g..lLarson. Tinslev fc Caldwell Il98d : 
iBalogh. Navarro fc Morris! l2000l : iMcCarthv et al.l l2008bF . 
The latter process removes the possibility for gas to cool 
and replenish the cold disc, leading to star formation fading 
away over a period of several Gyr and is therefore commonly 
called 'starvation' or 'strangulation'. This process is widely 
believed to be (at least partially) responsible for the the 
obser ved relation between galaxy colour an d environment 
(e.g.. lWeinmann et al.ll2006l : lFont et all2008l ). For X-ray ob- 
servations, this implies less massive hot gas haloes around 
galaxies in groups and clusters, and therefore lower X-ray 
luminosities, compared to galaxies of similar mass in the 
field. 

A different effect of the ICM on the hot gaseous 
coronae of galaxies has recently been proposed by 
iMulchaev fc Jeltemal (|2010l . MJ10). Motivated by observa- 
tions which appear to indicate a relative excess of X-ray 
emission by galaxies in groups and clusters, compared to 
galaxies of the same K-band luminosity in the field, these 
authors suggested that the very hot (Ticm ~ 10 K) ICM 
could exert pressure on cluster galaxies and their relatively 
cool gas halos with Thaio ~ 10 6 K. This 'confinement pres- 
sure' could prevent the outflows driven by supernovae, AGN 
or massive stars from leaving the galaxy and therefore poten- 
tially result in an increase in the hot gas density in galaxies 
orbiting within groups and clusters compared to those of 
similar (stellar) mass in the field. 

In this paper, we aim to investigate the rela- 
tive importance of ram pressure stripping and confine- 
ment pressure on group and cluster galaxies using the 
Galaxies-Inte rgalactic Medi um Interaction Calcu- 
lation (gimic: ICrain et al.l 120091 . C09), a set of five high- 
resolution cosmological hydrodynamical re-simulations of 
the formation and evolution of galaxies in a wide range of 
large-s cale environments ta ken from the Millennium Simu- 
lation jSpringel et al.l [20051 ). Analysing simulated galaxies, 



rather than real ones, brings the key advantage that the 
physical quantities responsible for both ram and confine- 
ment pressure, particularly the mass-weighted temperature 
and galaxy 3D velocity, are easily available. A second benefit 
is the possibility to trace the evolution of individual galaxies 
over time, which, as we demonstrate in Section [31 is key to 
understanding the potential effect of confinement pressure. 
Furthermore, a particular advantage of using GIMIC is that 
its field galaxies have already been shown to have properties 
in good agreement with observational data, such as scaling 
of X-ray luminosity with K-ba nd luminosity, star -formation 
rate and disc rotation velocity (|Crain et al.l2010ah and those 
of stellar spheroids around Milky Way mass disc galaxies 
|Font et al.ll201ll : lMcCarthv et alj|2012j ). Analysing a real- 
istic simulation will thus allow us to make meaningful con- 
clusions concerning galaxies in the real Universe. 

This paper is structured as follows: In Section 2 we de- 
scribe the cosmological simulations and our data extraction 
method, and present our results from both GIMIC and a sim- 
ple analytic model in Section 3. Our findings are summarised 
and discussed in Section 4. A flat ACDM cosmology with 
Hubble parameter h = Ho/XlOOkms^Mpc" 1 ) = 0.73, dark 
energy density parameter J7a = 0.75 (dark energy equation 
of state parameter w = — 1), and matter density parameter 
Qm = 0.25 is used throughout this paper. 



2 SIMULATIONS AND ANALYSIS 

2.1 Simulations and sample selection 

We extract galaxy groups and clusters from the Galaxies- 
Intergalactic Medi um In teraction Calculation suite 
of simulations (gimic: |C09| ). a set of five high-resolution (a 
baryon mass resolution of m gas ~ 1.16 x 10 7 /i _1 Mq with a 
gravitational softening that is 1 hT x kpc in physical space 
at 2 < 3 and is fixed in comoving space at higher red- 
shifts) re-simulations of nearly spherical regions of vary- 
ing mean density extrac ted from the Millennium Simula- 
tion (|Springel et alj|2005h . The simulations were carried out 
wi th the T r eePM -SPH code GADGET-3 (last described 
in ISpringell 120051') and include presc riptions for star for- 



mation ( Schave fc Dalla Vecchiall2008h. metal-dep endent ra- 
diative cooling (|Wiersma. Schave fc Smith] I2009T ). feedback 
and mass transport by Type la and Type II supernovae 

as well as stellar evolution 
The reader is re- 
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for a detailed description of the simulations. 

As we are interested in groups and clusters, we u se th e 
two highest density simulations, ' + 1<j' and '+2cr' (see lC09h . 
Groups and clusters of galaxies were identified at redshift 
z — using a standard Friends-of-Friends (FoF) algorithm 
with linking length b = 0.2 times the mean inter-particle sep- 
aration. We select all FoF groups with M200 > 10 13 M Q , 
where M200 is the mass within a spherical region of ra- 
dius V2oo, centered on the most-bound particle, in which the 
mean density is 200 times the critical density of the universe. 
There are 69 systems in total in the +lo and +2er gimic re- 
gions, with masses in the range 13.0 < log 10 (M2oo/A^e) < 
15.1. We explore below whether our results depend on total 
group/cluster mass, M^oo- 
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Within the simulated groups and clusters, bound sub- 
structures a r e ide ntified using the SUBFIND algorithm of 
iDolag et all (l2009h. that exte nds the standard implemen- 
tation of Springel et al.l (|200ll ) by including baryonic parti- 
cles in the identification of self-bound substructures. 'Galax- 
ies' are identified as self-gravitating substructures with to- 
tal stellar mass of Af* > 10 9 Mq (i.e., similar to the mass 
of galaxies typically identifiable in present observations of 
local groups and clusters). We exclude the central domi- 
nant galaxy (the 'BCG') from our analysis, as we are in- 
terested in the competition between ram pressure stripping 
and confinement pressure on orbiting galaxies. Also excluded 
are 'galaxy' identifications at r < 0.02 7-200 which we found 
to be associated with transient substructure in the BCG 
rather than being independent objects. In the 69 groups 
and clusters we have selected there are 1447 galaxies that 
meet our selection criteria, with stellar masses in the range 
9.0 < log 10 (M„/M Q ) < 12.3. 



2.2 Analysis 

For each of these galaxies we compute the ratio of confine- 
ment pressure to ram pressure at the position of the galaxy 
as follows. The ram pressure is given by P ram = pv 2 , where 
p is the density of the surrounding intracluster or -group 
medium and v the velocity of the galaxy relative to the ICM. 
The confinement pressure is simply the thermal pressure ex- 
erted by the hot ICM and is thus given by the ideal gas law 
-Pthermai = nk B T where n is the intracluster particle number 
density, k B Boltzmann's constant and T the temperature of 
the ICM surrounding the galaxy. Combining these two, we 
obtain 



pm p v 2 



(1) 



where p is the mean molecular weight of the ICM and m p 
the proton mass. We compute the galaxy velocities v relative 
to the mass- weighted average velocity of the group or cluster 
particles which implicitly assumes the ICM is at rest in the 
cluster centre of mass frame; for simplicity we assume a con- 
stant value of p = 0.58 throughout. For each group and clus- 
ter, we compute a spherically averaged radial temperature 
profile by binning the hot gas (T > 10 5 K) particles in bins 
of radial width A r = 0.1 r2oo- The ICM temperature at the 
position of each galaxy is then determined by linear interpo- 
lation. We have verified that using locally-determined ICM 
temperatures computed as the mass-weighted mean of all 
hot gas particles within a radial range 20 kpc < r < 50 kpc 
around the galaxy centre produces very similar results. 

Note that the ratio a between confinement and ram 
pressure in equation Q] depends only on the galaxy velocity 
v and the ICM temperature T, it does not depend on the 
density distribution of the ICM or the structure of the galaxy 
itself. The velocity, in turn, depends only on the depth of the 
potential well of the group/cluster, which is dominated by 
dark matter, and to fir st order this is also true of the temper - 
ature of the ICM (e.g.. lVoit et al.ll2002l : [Hansen et al.ll201lf ). 
Therefore, we expect our results to be insensitive to uncer- 
tain sub-grid processes such as star formation and feedback. 
Furthermore, since the dark matter mass distribution is ap- 
proximately self-similar, we also expect our results to be 



approximately independent of total group/cluster mass or 
redshift. We explicitly verify this below. 



3 RESULTS 

3.1 Analytic expectations 

Before proceeding to an analysis of the cosmological hydro- 
dynamic simulations, we can gain some insight by consider- 
ing a simple spherical analytic cluster in which the cluster 
galaxies orbit with the typical velocity v c i IC (r) = [GM tot (< 
r)/r] 1//2 and the hot gas is in hydrostatic equilibrium within 
the cluster potential well. 

Re-writing the equation for hydrostatic equilibrium, 



dP 

dr 



GM tot (< r)p 



dlnr - " clrc(rJ P' 



and using the ideal gas law 
P = ■ 



pm p 



-k B T 



(2) 
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with mean molecular weight p and proton mass m p , we ob- 
tain 



dlnP 2 ("n P 

= — Wcirc(r) 



1 

a 



(5) 



dlnr v ' k B T 

where the last equality is from equation ([T]) above, assuming 
galaxies move at the circular velocity v c irc(r). The typical 
ratio between thermal and ram pressure is therefore directly 
related to the slope of the logarithmic pressure profile of the 
host group or cluster. 

Furthermore, if the hot gas density distribution is as- 
sumed to trace that of the dark matter and both are de- 
scribed by a spherically symmetric power-law of the form 



-P 

p oc r 

then the mass enclosed within a radius r follows 



f 2 
M(< r) = / 47rr pdr 

Jo 

From equation ((2]), we then obtain 
dP _ GM(< r) 



oc r 



dr 1 

so that 

P oc r 

and therefore from equation (5) 
'dlnP 



-p oc 



2-2^ 



d In r 
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(6) 



(7) 



(8) 



(0) 



(10) 



For the case of an isothermal sphere with /3 = 2, the 
ratio of confinement pressure to ram pressure is a — 1/2. For 
a more realistic NFW profile, in which the effective exponent 
varies between /3 = 1 in the innermost regions and /3 = 3 
at large distances from the centre, a should vary as well: 
Equation (| 10p predicts asymptotic behaviour of a — ► 00 for 
r — ¥ 0, and a —¥ 0.25 for r —¥ 00, with a = 0.5 at the cluster 
scale radius (~ 0.20 — 0.25 r2oo)- 
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These analytic calculations therefore suggest that con- 
finement pressure will generally be subordinate to ram pres- 
sure except in the very inner regions of galaxy groups and 
clusters. However, there are several potentially important 
caveats to the above argument. First, the assumption that 
galaxies move at a velocity of exactly v c i IC (r) is clearly not 
correct and dispersion in velocity will cause scatter in the 
pressure ratio a which might lead to confinement pressure 
being important beyond the very centre. Secondly, devia- 
tions from spherical symmetry and hydrostatic and virial 
equilibrium may be relevant. Finally, the gas distribution 
will not follow that of the dark matter precisely and this 
will have an effect (albeit a small one) on the temperature 
of the ICM. 

For an as realistic as possible answer to the question of 
whether confinement or ram pressure stripping is the more 
important influence of the ICM on galaxies, we therefore 
need to take these complicating factors into account as well. 
This requires use of a detailed hydrodynamic simulation of 
groups and clusters, such as GIMIC. 
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3.2 Relative importance of confinement and ram 
pressure in GIMIC 

Fig. [T] shows the pressure ratio a obtained as described in 
Section [2] as a function of the cluster-centric radius r of each 
galaxy, normalised to the virial radius V2oo of the host. Each 
of the filled grey circles represents a simulated galaxy, while 
the black open diamonds and error bars give the median and 
25 th /75 th percentile of the distribution within bins of width 
A(r/r 200 ) = 0.2. 

For a clear majority of galaxies, the ratio is less than 
unity, and so ram pressure is dominant. But despite this 
general result, there are individual galaxies for whom con- 
finement pressure exceeds ram pressure (16 % of the whole 
sample at z = 0), in some cases by more than a factor of 
ten. 

For comparison, we construct somewhat more realis- 
tic analytic cluster models than considere d in Section 13.11 
In part icular, we use the methodology of iMcCarthv et alj 
(|2008al ) to generate clusters in hydrostatic equilibrium 
within an NFW potential well for several different central en- 
tropy levels: (yellow curves), 30 keV cm 2 (red curves) and 
100 keV cm 2 (blue curves). These roughly span the range 
of observed central en tropies of local galaxy clusters (e.g., 
ICavagnolo et a.lj|2009h - The solid colour curves correspond 
to galaxies orbiting with the typical velocity of v c i IC (r). The 
dotted colour curves correspond to case where galaxies or- 
bit with a velocity of Vci IC (r) + o"3d(»"), where (J3n(r) is the 
3D velocity dispersion profile of the NFW cluster. This case 
represents an effective lower limitQ on the ratio of the con- 
finement pressure to the ram pressure. 

There is remarkably good agreement in the typical ratio 
seen in the simulations with what is predicted by the ana- 
lytic models at all radii except for the very smallest. In both 

1 Note that since (r) is generally larger than v c i IC (r) we can- 
not compute an analogous upper bound curve. Physically, the up- 
per bound occurs at apocentre, when the galaxy's orbital velocity 
is at a minimum. For purely radial orbits, the ratio of confinement 
pressure to ram pressure goes to infinity at apocentre. 



Figure 1. Ratio a of confinement (thermal) pressure to ram pres- 
sure for our sample of galaxies from the GIMIC simulation, plotted 
versus the galaxies' cluster-centric distance in units of T200 (grey 
circles). The black diamonds and error bars show the resulting me- 
dians and 25 th /75 th percentiles in bins of width A(r/r20o) = 0.2 
For most galaxies — but not all — ram pressure dominates over 
confinement pressure. The solid yellow, red and blue curves show 
the expected relation from a simple analytic model of gas in hy- 
drostatic equilibrium inside an NFW halo as described in the text. 
The corresponding dotted curves indicate the expected lower lim- 
its. The dashed red curve corresponds to the expected relation 
from a simple analytic model of gas in hydrostatic equilibrium 
inside an NFW halo that takes into account the mass distribu- 
tion of the BCG (with the corresponding lower limit shown by 
the red dash-dot curve). The horizontal green dashed and dash- 
dot-dot-dot lines indicate a ratio a of 0.5 and 0.25, as expected 
from a power-law density profile with exponent -2 and -3 respec- 
tively. Overall, the galaxies from the GIMIC simulation follow the 
predictions of the analytic models quite well. 



the simulations and the analytic models, the pressure ratio 
is largely independent of position within the cluster, except 
for small radii (r/r200 < 0.2), where the predicted pressure 
ratio from our model increases sharply and tends to infin- 
ity as r — > 0, as expected (see Section 3.1). The galaxies in 
GIMIC, however, show only a small increase. 

A possible reason for this discrepancy at small radii 
is the presence of the BCG in GIMIC. Its mass acts to in- 
crease the slope of the potential well in the cluster centre 
and therefore, according to the logic in Section 13.11 to re- 
duce the pressure ratio a in the central region. To investigate 
the extent of this effect, we have modified the 30 keV cm 2 
model (red curve in Fig. [TJ to take into account the presence 
of the BCG of the massive cluster at the centre of the +2a 
simulation (i.e., we recompute the hydrostatic configuration 
of the gas using the sum of the dark matter NFW and BCG 
mass distributions). The pressure ratio resulting from this 
modified model is shown as a red dashed curve in Fig. [T] and 
is in much better agreement with the simulation results. In 
any case, only a small fraction of galaxies spend time there 
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Figure 2. Variation of confinement to ram pressure ratio with 
halo mass. The three curves show the distribution functions for 
three halo mass ranges as indicated in the figure. The distribution 
of confinement to ram pressure ratio does not depend on halo 
mass. 

and those that do venture into the very centre will likely be 
disrupted by tidal forces. 

A comparison of the three analytic model curves also 
confirms our expectation that the results are quite insensi- 
tive to the entropy (or density) distribution of the ICM. 

The extent of variation of the ratio between confine- 
ment and ram pressure with halo mass is shown in Fig. f2] 
No systematic variation is evident, which is consistent with 
our expectations based on the self-similarity of dark matter 
halos. 



3.3 Effect of confinement domination 

Despite confirming the overall dominance of ram pressure 
stripping in groups and clusters, the above results show that 
there is a small fraction of galaxies in which confinement 
pressure is expected to be more important than ram pres- 
sure. However, this does not automatically imply that these 
galaxies contain more hot gas than ram -pressu re dominated 
ones (or, more relevant in the context of lMJldl . field galaxies 
of the same stellar mass) — the relative importance of both 
may well have changed since the galaxy joined the group 
or cluster. iMcCarthv et al.l (|2008bh found that ram pressure 
stripping of hot gas halos takes place mostly during the first 
~ 2-3 Gyr after infall: At this point, the galaxy has typi- 
cally passed the first pericentre of its orbit and lost all hot 
gas that can be stripped. In order to enhance the hot gas 
content of a galaxy, confinement pressure must therefore be 
dominant within this period, or there will not be much gas 
left to be confined. To determine whether confinement pres- 
sure will indeed act to retain hot gas haloes in group and 
cluster galaxies, we therefore need to track the orbital his- 
tory of the galaxies from the time of infall to the present 
day. 

For this purpose, we define 'infall' as the first time 
a galaxy crosses the radius r = 2r2oo(z), where we have 
adopted this radius, rather than e.g. r2oo, because many 



galaxies return to cluster-centric distances V2oo < r < 2 r2oo 
after pericentric passage (see, e.g., Fig. [3} and therefore en- 
ter the region within r2oo more than once. To determine the 
infall time of the galaxies identified at redshift z — we 
trace their dark mat ter halos ba c k in t ime by a method sim- 
ilar to that used by iFont et all (|201lh . Firstly, we identify 
for each 'galaxy' subhalo at redshift z = the dark matter 
particles associated with it by SUBFIND, and then use their 
unique simulation IDs to find the subhalo to which each of 
these dark matter particles belonged to at redshift z = 2, 
discarding unbound particlefl The subhalo that contains 
most of these DM particles (excluding the BCG) is iden- 
tified as the main galaxy progenitor. At each intermediate 
redshift z, we then find the subhalo containing most of the 
DM particles present in both our original halo (at 2 = 0) 
and its progenitor at z = 2 to trace the galaxy forward in 
time. The same procedure is applied to the entire FoF group 
itself, starting with all DM particles with n < r2oo at z = 
to find the accurate r2oo(z) in each snapshot. 

The infall time tin is then found by linear interpolation 
between its position at the last snapshot with r ga i > 2 r2oo 
and the one immediately after this. Although this method 
cannot identify the infall time for galaxies already in the 
cluster at z = 2, this only affects a small fraction (6%) of 
the galaxies in orbit about the groups and clusters at z = 
and is insignificant for our conclusions below. 

The top and middle panels of Fig. [3] show, respectively, 
the ratio a between confinement and ram pressure and the 
cluster-centric distance normalised to the present-day virial 
radius of the host group/cluster, in both cases as a func- 
tion of time since infall. Both show a clear trend with time: 
the distance decreases for the first ~ 2 Gyr after infall as 
galaxies approach the first pericentre passage, and increases 
again until ~ 4 Gyr as the galaxies get closer to their first 
apocentre. Correspondingly, the ratio a is very low (i.e., ram 
pressure is strongly dominant) for the first ~ 2 Gyr with no 
single galaxy being confinement-dominated and a median 
value a ~ 0.1. Only around 4 Gyr after infall, at the time of 
first apocentric passage and the associated drop in galaxy 
velocity is confinement becoming increasingly more impor- 
tant. However, by this point in time many of the galaxies 
have already been completely stripped of t heir hot gas (bot- 
tom panel; see also lMcCarthv et al.ll2008bh . The most active 
period of ram pressure stripping occurs near first pericen- 
tric passage and since first pericentre necessarily precedes 
first apocentre there will be much less (if any) hot gas to be 
confined at first apocentre. 

A close inspection of the bottom panel of Fig. [3] shows 
that even at infall only ~ 45% of simulated galaxies with 
log 10 (M*/Mq) > 9.0 have hot gas atmospheres. This rel- 
atively low fraction is due to two effects: mass selection 
and the influence of the group/cluster beyond 2r2oo- In 
terms of mass selection, in the simulations a galaxy of stel- 
lar mass log 10 (M*/M©) = 9.0 has a total virial mass of 
log 10 (A^20o/A/e) ~ 11.0 ± 0.2, which lies approximately on 
the halo mass t hreshold required to sup port a hot gas at- 
mosphere (e.g.. iBirnboim fc Dekell l200St ) . Therefore, even 
amongst the GIMIC field galaxies of this mass, ~ 25% have 



2 In our adopted cosmology, a redshift of z = 2 corresponds to a 
lookback time of ~ 10 Gyr. 
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Figure 3. The ratio a = Pthcrmal/^ram (top panel), cluster- 
centric distance (middle panel), and fraction of galaxies with some 
hot gas still bound (bottom panel) as a function of time since in- 
fall into the group or cluster, see text for details. Each filled grey 
circle represents a simulated galaxy. The filled black squares with 
error bars represent the median and 25 th /75 th percentile of the 
distribution within each bin. Galaxies that have been in the clus- 
ter for less than ~ 2-3 Gyr (corresponding to the time between in- 
fall and first pericentric passage) arc all ram-prcssurc dominated. 
Pressure confinement becomes increasingly important near apoc- 
entre (roughly 4 and 8 Gyr after infall). The bottom panel shows 
the fraction of galaxies with some hot gas still bound, which de- 
creases rapidly between 2 and 4 Gyr after infall. The error bars 



1 .0 




0.0 I , , , , 

0.1 1.0 10.0 
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Figure 4. Fraction of galaxies with bound hot gas, plotted 
against the ratio a = P, thermal/-R'am ■ The solid and dashed curves 
represent the full galaxy sample and those with Af* > 10 10 Mq, 
respectively. While the latter are, in general, more gas-rich, both 
samples show a clear trend to less hot gas with increasing impor- 
tance of confinement pressure. 



no hot gas haloes, whereas these are found in virtually all 
more massive field galaxies (log 10 (M,/M Q ) > 10.0). The 
second and more imprtant factor is that ram pressure strip- 
ping is effective in the simulations at radii well exceeding 
2r2oo, the radius whose crossing we define as 'infall' in this 
study. As we will show in a forthcoming paper (Bahe et al., 
in prep.), the hot gas atmospheres of massive clusters induce 
stripping out as far as ~ 4>2oo from the cluster centre. As a 
result, roughly half of the originally gas-rich galaxies reach 
2r2oo without any remaining hot gas. Interestingly, there is 
mounting observational evidence that the effect of environ- 
ment on_g_ala : xieji(h2ej^nd^ several virial radii 
(e.g., lHaines et al"1l2009l ; iLu et alj|20l3 ). 

As further confirmation of the general ineffectiveness 
of confinement pressure, we show in Fig. U the fraction of 
galaxies containing any bound hot gas (T > 10 5 K) as a 
function of the ratio between confinement and ram pressure. 
When ram pressure is most strongly dominating, this is the 
case for ~ 20% of galaxies. With increasing importance of 
thermal confinement pressure, however, this fraction actu- 
ally decreases, reaching a level of only ~ 10% out of those 
galaxies with Ahermai > fram- This trend is even stronger 
for more massive galaxies (M, > 10 10 Mq): within these, 
hot gas is present in ~ 60% at P t hcrmai/-Pram ~ 0.1, but 
only in ~ 20% of confinement dominated cases. Therefore, 
confinement pressure is not only generally inefficient com- 
pared to ram pressure, but actually appears counter-effective 
in keeping hot gas within the halo in those cases where it is 
dominant. This is due simply to the fact that galaxies that 
are currently pressure confined have already been heavily 
stripped (as shown above they have gone through first peri- 
centre), whereas galaxies that are ram pressure dominated 
may not yet have been fully stripped if they have fallen in 
only recently. 
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Figure 5. Thermal X-ray fluxes (filled symbols) and upper lim- 
its (open symbols) as function of galaxy K-band luminosity from 
the literature, see text for details. Considering the full sample 
of galaxies shown here, low-Ljj field galaxies do not appear sys- 
tematically less X-ray luminous than their group/cluster counter- 
parts. 



4 DISCUSSION AND CONCLUSIONS 

We have investigated the relative importance of ram pres- 
sure stripping and thermal pressure confinement on galaxies 
in groups and clusters at redshift z — in the GIMIC simula- 
tion. A large sample of 69 groups and clusters with masses in 
the range 13.0 < log 10 (M2oo/A^©) < 15.1, containing over 
1000 galaxies in total, was analysed for this purpose. Our 
findings may be summarised as follows: 

• Thermal confinement pressure only dominates ram 
pressure in a small fraction (16%) of galaxies. In the ma- 
jority of cases, the action of the ICM is a removal of hot gas 
from group and cluster galaxies. 

• The ratio between confinement and ram pressure ob- 
tained from GIMIC agrees well with results from a simple 
analytic model of the ICM in which gas is in hydrostatic 
equilibrium within an NFW potential well and galaxies or- 
bit at the typical velocity of Vd TC (r). Outside the very central 
cluster region, the median value of the ratio is a ~ 0.3. 

• Increased confinement pressure does not lead to in- 
creased retention of hot gas. On the contrary, the fraction 
of galaxies containing any hot gas decreases with increasing 
importance of confinement pressure (from ~ 60% to ~ 20% 
for galaxies with M* > 10 10 M Q ). 

• The ineffectiveness of confinement pressure to retain 
hot gas is explained by the orbital history of the galaxies: 
upon falling into the cluster, they first experience the max- 
imum ram pressure stripping influence around pericentric 
passage, and only later — when they reach apocentre — a 
dominant influence of confinement pressure. By this point, 
there is often no hot gas left to be confined. 



We point out that the ratio a between confinement and 
ram pressure depends only on the galaxy velocity v and 
the ICM temperature T. Both of these quantities are set 
by the depth of the cluster potential well, which is domi- 
nated by dark matter. Therefore, we expect our results to 
be insensitive to uncertain sub-grid processes. Furthermore, 
since the dark matter mass distribution is approximately 
self-similar, our results are approximately independent of 
total host mass, which we explicitly verified in Section [3.21 

Our results are therefore seemin gly at odds with the 
observational evidence presented by iMJld . which appar- 
ently show that galaxies of fixed stellar mass in groups and 
clusters are more X-ray luminous than their field counter- 
parts. A potential explanation for thi s is sh own in Fig. [5] 
an expanded version of Fig. 1 from iMJlfj . which shows 
the Lk — Lx relation for fi eld and group/clus ter galaxies. 
The f ormer are taken from IMJld . as well as iDavid et al.l 
(l2006l. D06), while we use the c atalogues pu blished by 
iJeltema. Binder fc Mulchaevl (|200Sl . JMB08) and lSun et all 
(|2007l . S07) for group and cluster galaxies respectively. All 
these catalogues contain in part galaxies in which a thermal 
X-ray component was not detected, and for which therefore 
only upper limits on the thermal X-ray flux are available 
(shown by open symbols in Fig. [S] We see no significant dif- 
ference between the X-ray properties of the detected galaxies 
in the field, group and cluster samples. 

However, in Fig. [S] we show the detection fraction s of 
grou p and cluster glaxies from the samples of IJMBoH and 
|S07I . In both samples galaxies with higher Lk — and there- 
fore deeper potential wells — are more likely to be detected 
in X-rays. Galaxies of a given stellar mass (as indicated by 
their K-band luminosity), on the other hand, are more likely 
to be detected in groups than in clusters. Both of these 
trends are what one qualitatively expects from ram pres- 
sure stripping: high mass galaxies are able to retain more 
gas (which accounts for the increased detection fraction with 
increasing Lk) and galaxies of fixed mass are more effec- 
tively stripped in clusters than in groups, due primarily to 
the higher orbital velocities of galaxies in clusters. We note, 
however, that the comparison between the group and clus- 
ter galaxy samples is made complicated by the fact that the 
surface brightness of the background hot gas (the ICM) is 
generally higher in more massive systems. Defining the de- 
tection fraction to be the fraction of galaxies above a certain 
fixed surface brightness threshold (which could be that of the 
surface brightness of the most massive cluster in the sample) 
would be one way to rectify this problem. 

There is also an issue regarding the comparison of the 
group and cluster galaxy samples to the field galaxy sam- 
ple. In particular, the group and cluster galaxy samples are 
based on an optical selection, whereas the field galaxy sam- 
ples are drawn from a mixture of previous optical and X-ray 
catalogs. X-ray follow-up of an optically-selected field sam- 
ple is required to make a proper lik e-with-li ke c omp arison 
with the group and cluster results of IJMBfjil and ISOTl . 

Our theoretical results strongly suggest that ram pres- 
sure stripping is generally dominant over pressure con- 
finement and that strangulation should be effective in 
the group/cluster environment. This bodes well for semi- 
analytic models of galaxy formation that invoke strangu- 
lation to explain the environmental dependence of galaxy 
colours for galaxies of fixed stellar mass. In a future study, we 
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Figure 6. X-ray detection fractions for the observational samples 
shown in Fig. [5] Shown is the fraction of galaxies above a given K- 
band luminosity Lk with detected thermal X-ray emission. While 
galaxies with larger Lk are, in all environments, more likely to be 
X-ray detected, there is a clear trend to lower detection fractions 
in denser environments at fixed threshold Lk- This confirms that 
galaxies in groups and clusters are hot gas poor compared to the 
field population, as expected from ram pressure stripping. 



will examine the effect of strangulation on simulated galax- 
ies in GIMIC. 
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